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GLOBAL WEAK SOLUTIONS TO COMPRESSIBLE QUANTUM 
NAVIER-STOKES EQUATIONS WITH DAMPING 

ALEXIS F. VASSEUR AND CHENG YU 


Abstract. The global-in-time existence of weak solutions to the barotropic compress¬ 
ible quantum Navier-Stokes equations with damping is proved for large data in three 
dimensional space. The model consists of the compressible Navier-Stokes equations with 
degenerate viscosity, and a nonlinear third-order differential operator, with the quantum 
Bohm potential, and the damping terms. The global weak solutions to such system is 
shown by using the Faedo-Galerkin method and the compactness argument. This system 
is also a very important approximation to the compressible Navier-Stokes equations. It 
will help us to prove the existence of global weak solutions to the compressible Navier- 
Stokes equations with degenerate viscosity in three dimensional space. 


1. Introduction 


In this paper, we are interested in the existence of global weak solutions to the barotropic 
compressible quantum Navier-Stokes equations with damping terms 

pt + div(pu) = 0, 

( 1 - 1 ) 


(/9u)t + div(/9U (g) u) -h Vp'*' — div(pDu) = —rou — rip|uf u -|- KpV 
with initial data as follows 


V Vp ) 


p(0,x) = po(ic), (pu)(0,x) = mo(x) in (1.2) 

where p is density, 7 > 1, u(8)u is the matrix with components UjUj, Du = i (Vu -h Vu'^) 
is the sysmetic part of the velocity gradient, and is the d—dimensional torus, here 

d = 2 or 3 . The expression is called as Bohm potential which can be interpreted as 
a quantum potential. The quantum Navier-Stokes equations have a lot of applications, in 
particular, quantum semiconductors [ 5 ], weakly interacting Bose gases [ 9 ] and quantum 
trajectories of Bohmian mechanics m- Recently some dissipative quantum fluid models 
have been derived by Jiingel, see [TO]. The damping terms 

—rpu — rip|u|^u 

is motivated by the work of [T]. It allows us to recover the weak solutions to dLlD by 
passing to the limits from the suitable approximation. The most importance is that the 
existence of solutions for the system dni) studied in the current paper is crucial to show 
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the existence of weak solutions for the Navier Stokes equations with degenerate viscosity, 
see HU. Models with these drag terms are also common in the literature, see mmm- 

When ro = ri = K = 0 in (jl.ip . the system reduces to the compressible Navier-Stokes 
equations with degenerate viscosity fi{p) = vp. The existence of global weak solutions 
of such system has been a long standing open problem. In the case 7 = 2 in 2D, this 
corresponds to the shallow water equations, where p{t, x) stands for the height of the 
water at position x, and time t, and u(t, x) is the 2D velocity at the same position, and 
same time. For the constant viscosity case. Lions in [12] established the global existence 
of renormalized solutions for 7 > |, and Feireisl-Novotny-Petzeltova [U and Feireisl |7| 
extended the existence results to 7 > |, and even to Navier-Stokes-Fourier system. The 
first tool of handling the degenerate viscosity is due to Bresch, Desjardins and Lin, see 
[3|, where the authors deduced a new mathematical entropy to show the structure of the 
diffusion terms providing some regularity for the density. It was later extended for the case 
with an additional quadratic friction term rp|u|u, see Bresch-Desjardins m- Meanwhile, 
Mellet-Vasseur m deduced an estimate for proving the stability of smooth solutions for 
the compressible Navier-Stokes equations. 

When ro = ri = 0 in (II.ip . the system reduces to the so-called quantum Navier-Stokes 
equations. Up to our knowledge, there are no existence theorem of weak solutions for large 
data in any dimensional space. Compared to the degenerate compressible Navier-Stokes 
equations, we need to overcome the additional mathematical difficulty from the strongly 
nonlinear third- order differential operator. We have to mention that the Mellet-Vasseur 
type inequality does not hold for the quantum Navier-Stokes equations due to the quantum 
potential. Thus, there are short of the suitable a priori estimates for proving the weak 
stability. Jiingel m used the test function of the form pep to handle the convection term, 
thus he proved the existence of such a particular weak solution. In a very recent preprint, 
Gisclon-Violet [ 8 ] proved the existence of weak solutions to the quantum Navier-Stokes 
equations with singular pressure, where the authors adopt some arguments in m to make 
use of the cold pressure for compactness. Our methodology turns out to be very close to 
their paper. Actually, the authors of [ 8 ] mention that the existence can be obtained 
replacing the cold pressure by a drag force. 


The existence of weak solutions to dni), with the uniform bounds of Theorem m is 
crucial for the existence of weak solutions to the compressible Navier-Stokes equations 
with degenerate viscosity in 3D, see |14j . In that work, we started from the weak solutions 
to dni), that is, the main result of this current paper. Unfortunately, the version with 
the cold pressure proved in [ 8 ] , is not suitable for the result in [lU • On the approximation 
in HU, we need the terms ri/ 9 |upu and Kp{ ^^ ) for proving a key lemma. In particular, 
inequality m is crucial to prove the existence of weak solutions to the compressible 
Navier-Stokes equations in 3D. This estimate is from the term npS/ 


We can deduce the following energy inequality for smooth solutions of (II.ip 


-^(0 + [ [ /o|Bup dx dt-|-ro f ( |u|^ dx dt-|-ri f f p\u\^ dx dt < Eq, (1-3) 

Jo Jn Jo Jn Jo Jn 
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where 


and 


E{t) = E{p,\i){t) = Q/ 3 |u|^ + f dx, 

Eo = E{p,u){0) = Qpoluol^ + dx. 

However, we should point out that the above a priori estimate are not enough to show 
the stability of the solutions of (II.ip . in particular, for the compactness of p^. we have 
the following Bresch-Desjardins entropy (see mm) for providing more regularity of the 
density 


Qp|u +Vlnpp +-^ + - rologp^ dx + J j^\V p^ dx dt 


+ 


r f pivu 

Jo Jn 


— V^uP dxdt + K 


[ plV^logpl^ 
Jn 


dx dt 


(1.4) 


< ^/ 5 o|uoP + iVy/^p + - rolog_/9o^ dx + C, 

where C is bounded by the initial energy, log_ g = logmin(g', 1 ). 

Thus, the initial data should be given in such a way 

Pq^L^{Q), Po > 0, V^gL^(H), -log_po G 


mo E T^(H), mo = 0 if po = 0, 


|mo| 

Po 


E 


(1.5) 


We define the weak solution (p, u) to the initial value problem (II.1|) in the following 
sense: for any t E [0, T], 

• (| 1 - 2 I) holds in P'(H), 

• m and (11.41) hold for almost every t G [0,T], 

• m holds in P'((0,T) x H)) and the following is satisfied 

p E L°°(0,T;L^(H)), E L°°{0,T-L^n)), 

VypEL°°(0 ,T;L2 (H)), Vp? G L\0,T-L^Q)), 

E L2(0, T; L2(H)), ^pVu E L‘^{0, T; 4.2(0)), 
p3u E L^(0, T; L^(0)), u E L\0, T; 4.2(0)), 

^|V2 logp| E 4.2(0, r;L2(0)). 

The following is our main result. 

Theorem 1.1. If the initial data satisfy ()1.5p . there exists a weak solution (p, u) to (H])- 
(II.2p for any 7 > 1, any T > 0, in particular, the weak solution (p, u) satisfies energy 
inequality o, BD-entropy (11.411 and the following inequality: 


ll^/p||L2(0,^;i^2(^I)) + l|Vp4 \\L^{0,T-,L^{n)) < C, 


(1.6) 
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where C only depends on the initial data. Moreover, the weak solution {p, u) has the 
following properties 

puGC([0,r];4eafc(^^))> (VP)i€L2((0,r)xfi); (1.7) 

If we use {pk,Uk) to denote the weak solution for k> 0, then 

y/p^Uf,, ^ y/pu strongly in L‘^{{0,T) X II), as k —>■ 0, (1-8) 

where {p, u) in (HHI) is a weak solution to (fTTD-ra with K = 0. We remark the metric 

3 

space C{[0,T]] of function f : [0, T] —)• L'i'(O) which are continuous with respect 

to the weak topology. 

Remark 1.1. We will use (HH)- dni) in m to prove the weak solutions to (jl.ll) with 
I’D = = K = 0. In fact, inequality ()1.6I) is very crucial to prove a key lemma in M- 

Remark 1.2. The existence result contains the case with «; = 0, which can be obtained as 
the limit when k > 0 goes to 0 in (jl.ip . by standard compactness analysis. 

Remark 1.3. The weak formulation reads as 

/ pu ■ dx\lZo — / / pvnft dx dt — / / pu 0 u : Vif dx dt 

Jn Jo Jn Jo Jn 


f f p"^ divif dx dt — f [ pBw ■. Vif dx dt 

Jo Jn Jo Jn 

: —ro / / ufjdxdt — ri / / p\u\'^u'il> dx dt — 2 k / / Ay/p\/y/pfj dx dt 

Jo Jn Jo Jn Jo Jn 


(1.9) 


— K 


Ay/py/pdivtp dx dt. 

I 

for any test function ip. 


in 


2. Faedo-Galerkin approximation 


In this section, we construct the solutions to the approximation scheme by Faedo- 
Galerkin method. Motivated by the work of Feireisl-Novotny-Petzeltova [ 6 ] and Feireisl 
[7], we proceed similarly as in Jiingel [TT]. We introduce a finite dimensional space Xj\f = 
spanjei, 62 ,...., cat}, where G N, each e* be an orthonormal basic of which is also 

an orthogonal basis of We notice that u G G°([0, T]; Xw) is given by 

N 

u{t, ic) = ^ \i{t)ei{x), {t, x) G [0, T] x 0, 

i=\ 


for some functions Xi{t), and the norm of u in G®([0, T]; X^r) can be written as 

N 


u||cO([0,r];A„) 


sup ^|Ai(t)|. 

i6[0,T] 


And hence, u can be bounded in C^{[0,T];C’^{Q)) for any A; > 0, thus 


u|lcO([0,T];C'=(D)) ^ C'(A:)||u||c-0([o_t];L2(Q)). 
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For any given u E ^^([O, T]-,Xn), by the classical theory of parabolic equation, there exists 
a classical solution p{t,x) E (^^([O, T]; (7^(0)) to the following approximated system 

pt + div(pu) = eAp, p{0,x) = po{x) in (0,r) x 17 (2.1) 

with the initial data 


/?(0, x) = po{x) > p > 0, and po{x) E C'°°(17), (2.2) 

where > 0 is a constant. 

We should remark that this solution p{t, x) satisfies the following inequality 

inf po{x) exp- l|divu|Uoo(n) ds < < suppo(a;) exp^° l|divu||ioo(n) ds 

x&fi 


for all {t,x) in (0,T) x 17. By (|2.2p and (12.3p . there exists a constant Oq > 0 such that 

0 < 6*0 < p{t,x) < — for (7, x) E (0,T) x 17. (2-4) 

c'o 

Thus, we can introduce a linear continuous operator S : C'°([0, T]-,Xn) —)• (^^([O, T]; (7^(17)) 
by ^(u) = p, and 

||S'(ui) - 5'(u2)||c0([0,r];C'=(O)) < C{N,k)\\ui - U 2 ||cO([o,t];I, 2 (o)) (2.5) 

for any k > 1. 

The Faedo-Galerkin approximation for the weak formulation of the momentum balance 
is as follows 


/ pu{T)ipdx— / mQipdx + p / / Au-Aipdxdt— / / (pu ® u) : Vy? dx 

Jn Jn Jo Jn Jo Jn 

f [ 2pBu : Vp dx dt — f f p^Vipdxdt + p f f p~^^V(pdxdt 

Jo Jn Jo Jn Jo Jn 


dt 


+ 


+ £ 


Vp • Vvupdxdt = —ro 


>0 Jn 

fT 


// 

do Jn 

rT 


up dx dt — ri 


10 Jn 


p|upupdxd7 


2k / / A^/pV^/p^l^ dx dt — K / / Ay^y^divi/^dxdi + (5 / / pVA^ppdxdt, 

do Jn Jo Jn Jo Jn 

( 2 . 6 ) 


for any test function ip G Xj\f. The extra terms pVp”^*^ and dpVA®p are necessary to keep 
the density bounded, and bounded away from zero for all time. This enables us to take 
^ as a test function to derive the Bresch-Desjardins entropy. 

To solve (j2.6p , we follow the same arguments as in [U [71 [H] and introduce the following 
operators, given the density function p(7, x) E L^(17) with p > p > 0, here we choose 
p = 6 q. We define 

2R[p(7), ■] : Xn ^ X]^, < T)l[p]u, w >= pu ■ w dx, 

Jn 

We can show that 911[p] is invertible 

ll^~H/5)llL(A^,A:iv) ^ 


for u, re E X^. 
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where L{Xy, Xj\f) is the set of all bounded linear mappings from X’^ to Xj\f. It is Lipschitz 
continuous in the following sense 

ll^~^(pi) “ ^~^{P2)\\l{X^,Xn) - “ P2||l1(o) (2.7) 

for any pi and p 2 from the following set 

= {p £ inf p > u > 0.} 

For more details, we refer the readers to ISIZlttl]. 

We are looking for u„ G C'([0,T];X„) solution of the following nonlinear integral equa¬ 
tion 


UN{t)=Tl ^(S'(uAr))(t) 5JI[po](uo)+ / ^{S{uN),UN){s)ds) , 


( 2 . 8 ) 


where 


frt(S'(uAr), uat) = —div(pu7v <8) ujv) + div(pDujv) + — eVp ■ Vu^r + pVp — Vp' 


- roUAT - rip|uAr|^UAr npV 


Ayp 

Vp 


+ 6pVA^p, 


p = ^(uAr). 

Thanks to (|2.5p and (j2.7p . we can apply a fixed point argument to solve the nonlinear 
equation ()2.8p on a short time interval [0,T*] for T* < T, in the space C'‘^([0 ,T*];Xat). 
Thus, there exists a local-in-time solution (pAr,ujv) to ()2.1I) . (|2.8p . Observe that — 
norm and norm are equivalent on X^r. 

Differentiating (12.6p with respect to time t and taking (p = u^v, we have the following 
energy balance 

-E{pN,ui\f) + p f |AuAr|^dx-|- f pAr|Bu7v|^ dx + eJ f |A®pAr|^dx 

' Jq Jq Jq 


dt 
-|- £ 


^ ,0 

2 2 

n\ 


dx + ep / |Vp^fdx-|-ro / \u]\fpdx + ri / pjvlujvrda; 


/ |Vp 
Jq 

£ / PAT Ilog PAT I = 0, 

Jq 


(2.9) 


on [0,T*], where 

E{pN,viN) = ^ + :^ + ° f + ^IVAVaP) dx. 


and 


Eq{pn,\in) = Qpoluop + ^Po^° -k ^IVAVoP^ dx. 


Here we used the identity 


2pArV(^^^) = div (pArV^(log pat)) 


to yield 


f Ay/P^ 
IQ ^/^ 


I 


ApM dx = - pnV log Pat • V 


dx=i f pArlV^logpATpdx. 

^/PN ) ^ Jq 
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Energy equality (|2.9p yields 


fT* 

/ ||AuAr||^2 dt < Eo(pn,U7v) < CX). 

Jo 

here exists a const an 
0 < 6*0 < PN{t,x) < 


( 2 . 10 ) 

( 2 . 11 ) 


Due to dimXTv < oo and (12.31) . there exists a constant 6*0 > 0 such that 

1 

To 

for all t E (0, T*). However, this Oq depends on N and it is the same to 9o in (12.41) . Energy 
equality ( 12 .9h gives us 

fT* 

PN\uN\‘^dx < Eo{pn,un) < 00 , 


sup / 
te{o,T*)Jo 


and 


rT* 

/ pAr|DuArp dx dt < Eo{pN, Uat) 
40 


< OO, 


which, together with (I2.10h . (12.lib , implies 


sup (lluArllioo + IIVuatIIloo + ||AuAr||L“) < C{Eo{pn,'OLn),N), 
oe(o,T*) 


( 2 . 12 ) 


where we used a fact that the equivalence of and L°° on X]\f. By (12.Sp . (|2.7I) . (|2.1ip 
and (|2.12p . repeating our above arguments many times, we can extend T* to T. Thus 
there exists a solution {pN,U]\f) to (|2.1I) . (|2.8p for any T > 0. 

Here we need to state the following lemma due to Jiingel m- 

Lemma 2.1. For any smooth positive function p{x), we have 


and 


p\X‘^ log p\^ dx > - / 


/?|V^ logpl^ dx > - / I'Vp'il^dx. 


Proof. The above inequality of Lemma is firstly proved by Jiingel m- Here, we give a 
quick proof. 

We notice 


^■V‘^logy/p = y/p- = V^^/p - ® 


Vp 


Vp 


(2.13) 


thus 


f pIV^ log dx = f |V^y^|^dx+ f |2V/?4|^dx —2 f V^^/p 

J ^ J ^ J ^ "/Q 

= A + B-I, 

Eor I, we control it as follows 

f^Vp 


Tp ’ 


1 = 2 

= -2 


V Vp 


VyJp dx 


Jn 

Jn 


In y/P 


Vp 


dx. 
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Hence: 


21 = -2 


f [Z^^Alog Viodx < 2 VWD, 

Jfl y/P 


Vp 

where D = p\V‘^ log dx, and hence 


A + B = D + I<{1 + 6)D + -B, 

O 


-A + -B < D. 
7 8“ 


and thus, 

So we proved this lemma. □ 

By (122]), we have 

E{PN,11 n) < Eo{p]y,UN), 

this gives us 

l|PA||l,°°(0,T;H9(n)) < C{Eo{pn,Un),S), 

this, together with ( 12 . 111 ) . gives us that the density p{t,x) is a positive smooth function 
for all {t,x). We also notice that 

rT 


KE 


/ / /OwlV^logpATpdxdt < £^o(/OA,UAr) 

Jo Jn 


< oo. 


By Lemma EH we have the following uniform estimate: 


(2.14) 


(k£)2 ||v^^||l 2(0,T;H2(O)) + (Ke)4 ||Vp^||2,4(o,r;L4(n)) < C, 

where the constant C > 0 is independent of N. 

To conclude this part, we have the following lemma on the approximate solutions 
{pN, uat): 

Proposition 2.1. Let {p]y,Uf^) be the solution of (|2.ip . (12.8p on (0,T) x Q constructed 
above, then we have the following energy inequality 

sup ujv) + ^ f f I AutvP da: (it + f f pn\^'^n\‘^ dx dt + eS ( f \A^p^l"^ dx dt 

te(o,T) Jo Jn Jo Jn Jo Jn 

+ £ [ [ \Vdx dt + ep f f I cix cii + ro f f dx dt 

Jo Jn Jo Jn Jo Jn 

PAr|u7v|^ da: (ii + Kff / / log pwl^ da: di < LioPPA, ujv), 

Jo Jn 


>0 Jn 


(2.15) 


where 


E{pn,um) = ^ QpaIujvP + :^ + ° f iVv^P + ^IVAVaP) dx. 

Moreover, we have the following uniform estimate: 

1 _ ]_ i 

{ k £)2 ||v^^||l2(o,T;A 2(0)) + (Ke)4 || Vp^ || A4(o,r;L4(n)) < C, 


(2.16) 
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where the constant C > 0 is independent of N. 

In particular, we have the following estimates 

^/PnUn € L°°(0, T; L^(r2)), y/pffDuN G -^^^((0, T) x H), y/JlAuisf G L^((0, T) x Q), (2.17) 

VldA^PN G L2((0,r) X n),V6pN G L°°{0,T- {!})), ^^ e L°°{0,T;H\n)), 

(2.18) 

VpIgL\{0,T)xQ), p-^ GL°°{0,T;L^\n)),^Vp],^ GL^mT)xn), (2.19) 
U 7 V G l2((0,T) X 12), G L^((0,r) X 12). (2.20) 

Based on above estimates, we have the following estimates uniform on N: 


Lemma 2.2. The following estimates hold for any fixed positive constants e, p, p and 6: 
ll(\/^)t||L2((o,r)xO) + ll\/^llL2(0,T;iL2(O)) < K, 

||(pAf)t||L 2 ((o,T)xn) + ll7'Af||L2(o,T;iLiO(0)) < K, 

\\{pNUN)t\\L^{0,T-,H-^{n)) + ll7’AfUAr||L2((o,r)xn) < 

V(/?ArUAr) is uniformly bounded in L^(0, T; Ls (12)) + L^(0, T; L 2 (12)). 

WpJ 

where K is independent of N, depends on e, p, p and 5. 

Proof. By ()2.16p . we have 

ll\/^llL2(o,T;iL2(r2)) < C. 


^^"L§(( 0 ,T)xr 2 ) “ 


"l§((0,T)xO) “ 


We notice that 


which gives us 


{PN)t = -pNdivUN - VpN ■ UAT 
i i 1 

= -{dVpff){pffUN)ipff) - ^/pff^/pffdlVUN, 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 


ll(PAf)4||L2((0,T)xn) < 4||V/9^||i4((o,r)xn)l|P^UAr||i4((o,T)xn)ll7’^llL°°((0,T)xO) 

+ ll^/^llL°°(( 0 ,T)x^ 2 )ll^/^VuJv||^ 2 ((o,T)x^I)) 

thanks to ()2.17p - (l2.20p and Sobelov inequality. 

Meanwhile, we have 

2 (v^)t = -^/pNdivu.N - 2Vy/pff ■ UN 

1 1 

= -y/pNdivUN - SVpffPf^UN, 

which yields (,ypN)t is bounded in L^((0,T) x 12). 

Here we claim that (/OatUjv)* is bounded in L^(0, T; i7“®(12)). By 

(pNUN)t = -div(pArUAr (g) Uat) - Vp]y + pVpfi^^ + pA^UN + div(/9ArDuAr) - PoUat 
- ri/37v|uArpUAr + eVPN ■ Vuat + kpnV 


+ dpNVA^pN, 


y/PN J 


we can show the claim by the above estimates. 
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And 

3 1 

l|PAfUAr||j^2((o,r)xQ) < l|PAr||L°°(0,r;L4(Q))||/J^UAr||i4((o,T)xr2) < AT, 
where we used Sobelov inequality and (I2.16p . Thus we have (I2.22I1 . 

We calculate 

1 1 

V{PNUn) = Vy/pNPnUNP^ + ^/pN^/m'^UN, 

it allows us to have (I2.23p . For any given e > 0, we have 

2 

II^PnIIr2((o,t)xD) < A', 

which gives us 

\\PN\\L^(0,T-,L^{n)) < A'. 

Notice 

pleL^{o,T-L\n)), 

we apply Holder inequality to have 


11 /^: 


N\ 


L3((0,T)xn) 


2 3 


Similarly, we can show (I2.25h . 

Applying Aubin-Lions Lemma and Lemma 12.21 we conclude 

Pn ^ P strongly in L^(0, T; weakly in L^(0, T; 

^J~Pn ^ ^fp strongly in L^(0,T; weakly in L^(0, T; 

and 

Ptvujv pu strongly in L^((0,T) x H). 

We notice that u^r E L^((0,T) x H), thus, 

U]sr —)■ u weakly in L^((0,T) x H). 

Thus, we can pass into the limits for term pn^n ® as follows 

Patutv <8) utv — >■ pu (g) u 

in the distribution sense. 

Here we state the following lemma on the convergence of pAr|uArpUAr. 
Lemma 2.3. When N oo, we have 

PAr|u 7 vpU 7 v p|u|^u strongly in L^(0,T; L^(H)). 
Proof. Fatou’s lemma yields 

/ p\u\'^ dx < / liminf patIuatI"^ dx < liminf / p7v|u7v|^dx, 

Jq, Jq, Jq, 

and hence p|u|^ is in L^(0, T; L^(H)). 

By (|2.26l) and (|2.27l) . we have, up to a subsequence, such that 

Pn ^ p(t,x) a.e. 


□ 


(2.26) 


(2.27) 


and 


Patutv —>• pu a.e. 
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Thus, for almost every {t,x) such that when pN(t,x) / 0, we have 

Pn^n 


UN = 


U. 


Pn 


For almost every {t,x) such that pN{t,x) = 0, then 

PN\uNfuNX\uN\<M < M^Pn = 0 = p\ufuX\u\<M- 

Hence, /5Ar|ujvpU7vX|ujv|<AL converges to p|upux|u|<^ almost everywhere for {t,x). Mean¬ 
while, PAr|ujvpU7vX|ujv|<M is uniformly bounded in L°°(0, T; L^(H)) thanks to (I2.18I1 . 
The dominated convergence theorem gives us 

PN\-iiNfuNX\uN\<M^pH‘^'^X\u\<M strongly in L^(0,r;L^(H)). (2.28) 

For any M > 0, we have 

/ / |/3Ar|uArpujv —/o|u|^u| da;dt 

Jo Jn 


< 


+ 2 


< 


/ / |/3Ar|uArpUArX|u^|<M -/J|uPuX|u|<m| dxdt 

Jo Jn 

/ / PN\uN\^X\uN\>Mdxdt + 2 / / p\ufx\u\>Mdxdt 

Jo Jn Jo Jn 

[ [ \pn\un\‘^unX\un\<m - pH‘^uX\u\<m\ dxdt 

Jo Jo. 


10 JQ 

Thanks to ()2.28p . we have 


+ ^ / / /aAr|uAr|'^da;dt -h 




dx dt. 


C 


lim sup ||/97v|u7v|^UAr -p|u|^u||il(o,r;Ll{n)) < TT 


for fixed C > 0 and all M >0. Letting M —oo, we have 

PAr|uAr|^UAr —>• p|upu strongly in L^(0,T; L^(H)). 


(2.29) 


□ 


By ()2.24l) and pj^ converges almost everywhere to p^, we have 

p]^ —)• p^ strongly in L^((0,T) x H). 

Meanwhile, we have to mention the following Sobolev inequality, see mm, 

\\p ^IIl°°(o) < ^(1 + lldllfi''+2(o))^(i-+ Up ^IIl3(o))^, 

for /c > |. Thus the estimates on density from (I2.17h - (l2.19h enable us to use the above 
Sobolev inequality to have 

IIpIIl'=°((o,t)xO) > C{6,p) > 0, a. e. in (0,T) x Q. (2.30) 

This enables us to have converges almost everywhere to p~^^. Thanks to (|2.25p . we 
have 


^ strongly in L^((0,T) x H). 
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By the above compactness, we are ready to pass into the limits as —)• oo in the 

approximation system (EH), ESI)- Thus, we have shown that (p, u) solves 

pt + div(pu) = eAp pointwise in (0, T) x n, 

and for any test function cp such that the following integral hold 


In 


pu{T)ipdx — / rriQipdx + p 

Jn 


0 Jn 


Au • Aif dx dt — 


/ / {pwS> u) : V(f dx dt 

Jo Jn 


+ 

+ £ 


2/9Du :V(pdxdt — 


fl 


p'^Vip dx dt + rj 


fl 


p ^^Vpdxdt 


V/9 • Vuipdxdt = —ro 


u(/? dx dt — ri 


lo Jn 

rT 


lo Jn 


p\u\‘^u(pdx dt 


0 Jn 

fl 

Jo Jn 

— 2k, f f A^/pV ^/p^l^ dx dt — K [ f A^/p^/pdivi/j dx dt + 6 f f pVA^ppdxdt. 
Jo Jn Jo Jn Jo Jn 

(2.31) 

Thanks to the weak lower semicontinuity of convex functions, we can pass into the limits 
in the energy inequality (I2.15p . by the strong convergence of the density and velocity, we 
have the following energy inequality in the sense of distributions on (0, T) 


sup E{p,u) + p f [ {Aul"^ dxdt+ ( 
E(o,T) Jo Jn Jo 

f 


uf dx dt + e5 


lo Jn 


A^pp dxdt 


+ £ 


+ ri 


2 p dx dt + £7? 


|V/? “fdxdt + ro 


luP dx dt 


lo Jn 

rT 


lo Jn 


p\u\^dxdt + Ke f f /?|V^logpp dxdt 

Jo Jn 


Jo Jn 
< Eo, 


(2.32) 


where 

In ^ 7^ ^ ^ f dx. 

Thus, we have the following Lemma on the existence of weak solutions at this level 
approximation system. 

Proposition 2.2. There exists a weak solution {p, u) to the following system 
Pt + div(/9u) = eAp, 

{pu)t + div(/9U (g) u) + — div(/9Bu) — /lA^u + eVp ■ Vu 

= -rou - ri/?|upu + KpV + dpVA^p, 

with suitable initial data, for any T > 0. In particular, the weak solutions (p, u) satisfies 
the energy inequality ()2.32l) and (I2.30p . 
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3. Bresch-Desjardins Entropy and vanishing limits 

The goal of this section is to deduce the Bresch-Desjardins Entropy for the approxima¬ 
tion system in Proposition 12.21 and to rely on it to pass into the limits as e, /i, rj, 6 go to 
zero. By ()2.18p and ()2.30p . we have 

p{t,x)>C{6,rj)>0 and p G L2(0,r;Fi°(D)) nL°°(0,r;F®(D)). (3.1) 


3.1. BD entropy. Thanks to we can use ip = V (log p) to test the momentum 

equation to derive the Bresch-Desjardins entropy. Thus, we have 

Lemma 3.1. 

SI (pl“ + T'" ^ ’'/n 

-|- f \Vpi\dx + 6£ f \A^p\‘^dx + 26 f \A^p\^dx + \- f p|Vu — V'^u|^ dx 

Jn Jn Jn 2 Jq 


+ p j \A\i\^ dx + K / p|V^ logpp dx-|-e 


f |Ap|- 
In P 


dx 


f f iVlogpP 

= e / Vp-Vu-Vlogpdx-|-e / Ap -^- dx — e 


— p An • VAlogpdx — ri / |u|^uVpdx —ro 

J ^ "/fi ■ 

= i?i -|- i?2 + d?3 -|- i?g. 


/H 2 

| 2 , 


/ div(pu)- 
Jn I 


■Apdx 


r u Vp 

'n P 


dx 


We can follow the same way as in |16] to deduce the above equality, and control terms 
Ri for i = 1,2,3,4, and they approach to zero as e —)• 0 or p —)■ 0. We estimate R^ as 
follows 

|d? 5 | < C f p|up|Vu| dx < C f p|u|^dx-|— f p|Vupdx, 

Jn Jn 8 Jq 

and for Rq we have 

Rfi = ro f ddivu —dx = xq f (log p)t dx - erg f — dx. 

Jn P Jn Jn P 

since p is uniformly bounded in L°°(0, T; L'''(D)), we have 

ro / log_,_ pdx < C, where log_,_ g = logmax(p, 1). 

Jn 

Thus, we need to assume that —xq J^log_podx is uniformly bounded in L^(D). Also we 
can control 

and it goes to zero as e ^ 0 . 

Thus, we have the following inequality 


£Xo 




In P 
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- rolog/ 1 ) dx + r? 


'o in 


|V/3 ^fdxdi 


/7 


+ 


< 


IV/?21 dxdt + (5e f f |A^/?|^ dxdt + 25 f f \A^p\‘^dxdt 
Jo in Jo in 

- [ [ /o|Vu — V^up dxdt +/X f f |Aupdx + K f p| logpp dx di 

2 io in io in in 

Vd2i+ e||/3||^^2(Q) 11/9"^ 11^00(0)+ C / //9|u|'^dxdt + ^ / //9|Vupdxdt+ 

Jo in o io in 

[ fJ/5o|uo + ^P + ^|VVoP + ^|VV^P + -^ + ^^-rolog.po) dx, 

in V2 Po 2 2 7 - 1 10 ) 

(3.2) 

where p|u|^ dx dt is bounded by the initial energy, and | p| Vup dx dt can be 

controlled by 


[ [ p|Vu-V'^u| 

io in 


dx dt, 


and 


/7 


up dx dt. 


In-deed, it can be controlled by 

rT 


J j p|Vupdxdt< J (^po\\io\'^ +-^^ + \V - rolog_ po'^ dx + 2Eo. 

Thus, ()3.2p gives us 


ip|u + —\‘^ + + 7vVpP -b - rologp ) dx -b ?? 


+ I I IV /921 dx dt-b de / / | A^pp dx dt-b 2d 


7-1 10 
T 


/o in 


|V/3 ^pdxdt 


A^yOp dx dt 


Io in 


'0 in 


/o in 


dxdt + p / |Aupdx-bxt / /9|V^ logpp dx dt 


/o in 


+l r [ pivu-v^ui 

2 = io in 

< 2 /" (^Po|uo + ^P + ^|VVoP + 7v^p-b-^-b^^-rolog_/9o) dx 

in V2 po 2 2 7 - 1 10 J 

4 

+ ^ -Ri + e||p||_H 2 (n) Up ^llL°°(n) + 2 T'o- 


Z =1 


Thus, we infer the following estimate from the Bresch-Desjardins entropy 

7 

K 


(3.3) 


; f f /9|V^ logpp dxdt < C, 
io in 


where C* is independent on e, p, p, d. 

Applying Lemma l2.11 we have the following uniform estimate: 

ll^/p||L2(0,^;i^2(n)) + l|Vp4 ||L4(0,r;L4(Q)) < C, 
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where the constant C > 0 is independent on e, t], jj,, 5. 


3.2. Passing to the limits as e, — >■ 0. We use to denote the solutions at 

this level of approximation. It is easy to find that has the following uniform 

estimates 

G L‘^{n)), e i^((o,r) x n),^Aue,^ e x n), 

(3.4) 


v/iAVe,M G L^mr) X G L^{0,T-,H^{n)),^^GL^i0,T-,H\n)) 


p-i G L°°(0,r;L'°(II)),v^Vp-5 G L\{0,T) x II), 

G L‘^{{o,T) X G L^((o,r) X n). 

By the Bresch-Desjardins entropy, we also have the following additional estimates 


G L°°(0, T; L‘^{n)), V6A^p,,^ G L^O, T; L\n)), 


and 


Vpl^ G L2((0,r) X n), 0?Vp-5 G L2((0,r) X n)). 

Also we have the following uniform estimate 

1 1 1 

\\^yPe,n\\L^{0,T-,H^{n)) + 1'^'^ l|VP£,At||L4(o,r;L4(n)) < C, 

where the constant C > 0 is independent of e, p, p, 5. 

By Lemma EH one deduces 


(3.5) 

(3.6) 

(3.7) 


(3.8) 

(3.9) 

(3.10) 




dx dt < C, 


which together with (13.411 . yields 

[ [ Pe^fj.lVUe^f.l'^ dxdt < C, 

Jo Jn 

where the constant C > 0 is independent of e, p, p, 6. Based on above estimates, we have 
the following estimates uniform in e: 


(3.11) 


Lemma 3.2. The following estimates holds: 

||(•v/pH^)^llL2(0,T;L2(f^)) + \\y/f^\\L^{0,T-,H^{n)) < 

||(P£,At)t||^2(0T-L§(O)) + I|7 ’£,mIIl2(0,T;/I10(O)) < 

\\iPe,^lUe,^l)t\\L^{0,T■,H-^{n)) + \\Pe,^J,^e,^l\\L^{{0,T)xn) < 
V(p£^^Ue^^) is uniformly bounded m L^(0, T; Ls (II)) + L^(0, T; L 2 (H)). 

Wo'^ II 5 <K 

\\Pe,u lli§((o,r)xn) — 


(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 


where K is independent of e, p. 
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Proof. By (I3.4I) - (I3.11I) . following the same way as in the proof of Lemma 12.21 we can 
prove the above estimates. 

Applying Aubin-Lions Lemma and Lemma 13.21 we conclude 

Pe,r, ^ P strongly in (7(0, T;Lf®(n)), weakly in L^(0, T; 77^°(fl)), (3.17) 

^/pffl^^/p strongly in L^(0,T;iL^(O)), weakly in L^(0, T; iL^(n)) 

and 

Pe,^J.^e,^l /ou strongly in L^((0,r) x 12). (3.18) 

We notice that G L^((0,T) x 12), thus, 

—)• u weakly in L^((0,T) x 12). 

Thus, we can pass into the limits for term ® as follows 

Pe,r.^e,r. ® pu (g) U 

in the distribution sense. 

We can show 


Pe,fi\ue,^,\‘^Ue^mu ^ 2 >|upu strongly in L^(0,r;L^(12)) 
as the same to Lemma 12.31 

Here we state the following lemma on the strong convergence of y/p^Un, which will be 
used later again. The proof is essential same to m- 
1 

Lemma 3.3. IfpnUn is bounded in L^(0, T; L'^(12)), p„ almost everywhere converges to p, 
Priori almost everywhere converges to pu, then 

y/pfiU-n y/p^ strongly in L^(0, T; L^(12)). 


Proof. Fatou’s lemma yields 

/ p\u\^ dx < / lim inf pn,|u„|^ dx < lim inf / p^.lu^l'^dx, 

Jo. Jo. Jn 

and hence pluj"^ is in L^(0,T; L^(12)). 

For almost every {t,x) such that when pn(t,x) / 0, we have 


u„ = 


Pn^n 

Pn 


U. 


For almost every {t,x) such that pn{t,x) = 0, then 


\/Pn'^nX\un\<M — — 0 — \fp'^X\\i\<M ■ 

Hence, yffhLXinX\\in\<M converges to ^/pvLX\u\<M almost everywhere for (2, x). Meanwhile, 
v^WX|u„|<M is uniformly bounded in L°°{0,T-, L^{Q)). 

The dominated convergence theorem gives us 

\/^WX|u„|<M VpuX\u\<M strongly in L‘^{0,T-, L‘^{n)). 


(3.19) 
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For any M > 0, we have 

rT r 


10 JQ 


< 


+ 2 


< 


\y/f^Un - y/pu\‘^ dx dt 

[ [ - Vp^X\u\<m\‘^ dx dt 

Jo Jn 

/ / \VK^X\u„\>Mf dx dt + 2 / / \^ux\u\>Mf dx dt 

Jo Jn Jo Jn 

I I P'nJ^nX\\in\<M ~ \/P'^X|u|<m| dx dt 
Jo Jn 


(3.20) 


+ 


M 2 


'0 Jn 


Pn\'a.n\‘^ dx dt + 


M 2 


/ p|u| 

Jn 


^ dx dt. 


Thanks to ()3.19p . we have 


lim sup \\y/p^Un — y/])u\\L2i^Qx-,L^{n)) < 
e,fi—>-0 


M 


for fixed C > 0 and all M >0. Letting M ^ oo, we have 

V^u„strongly in L^(0,r;L^(n)). 

Applying Lemma 13.31 with (|3.17l) , (13.181) and 

rT 

dx dt < C < oo, 


□ 


/ / Pe,^l\^e,^l\ 

Jo Jn 


we have 

yfp^ strongly in L^(0, T; L^(n)). 

By (|3.15p and converges almost everywhere to p^, we have 

p2,ti / strongly in L^((0,r) x fi). 

Thanks to we have p^j^^ converges almost everywhere to p Thus, with p3.16p . 

we obtain 

Pe}^ strongly in L^((0,T) x fi). 

By previous estimates we can extract subsequences, such that 

P£,ii 0 strongly in L^((0,T) x 12), 

and 

-)■ 0 strongly in L^((0,T) x 12). 

For the convergence of term for any test function p € L^(0,T; M^(12)), we have 

rT 


/ / dx dt 

Jo Jn 


^ VP‘\\VpJ^''^s,IJ-\\L'^{0,T;L^{n))\\J^‘f\\L'^(0,T;L^{n)) 0 


as /U —>■ 0, thanks to p3.4p . 

Due to weak lower semicontinuity of convex functions we can pass into the limits in 
energy inequality (I2.32p . we have the following Lemma. 
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Lemma 3.4. 


+ 


/7' 


updxdt + ro / / \u\‘^dxdt + ri / / /?|u|' 

./o in io in 


dx dt 


(3.21) 


Passing to the limits in (13.3h as e —)• 0 and /i —^ 0, we have the following BD entropy. 

Lemma 3.5. 


In 


^p|u + ^1^ + -rologp ) dx+ 


' r [ ivp-'p 

io in 


dx dt + 


) 21 dx dt + 2(5 

rT 


[ |a7P 

in 


dx dt 


[ /|Vp^ 

io in 

+ 1: [ [ p|Vu — V^up dx dt + K f j p| logp|^ dx dt 

2io in io in 

[ f^Po|uo + ^P + ^|VVoP + ^|Vv/^|^ + -^ + ^^-rolog_po) dx 
in V2 Po 2 2 '' 7-1 10 ) 

(3.22) 


< 2 
+ ‘IEq. 


Thus, letting e ^ 0 and ^ 0, we have shown that the following existence on the 

approximation system. 

Proposition 3.1. There exists the weak solutions (p, u) to the following system 
pt + div(pu) = 0, 

(pu)t + div(pu (g) u) + Vp'^ — pVp“^° — div(pDu) 

= -xqu - rip|upu + KpV + (5pVA®p, 

with suitable initial data, for any T > 0. In particular, the weak solutions (p, u) satisfies 
the BD entropy (j3.22l) and the energy inequality (13.210 . 


3.3. Pass to limits as p, (5 —?> 0. At this level, the weak solutions (p, u) satisfies the BD 
entropy (j3.22j) and the energy inequality (I3.21|l . thus we have the following regularities: 

G L°°(0, T; L^{D)), G T) x D), (3.23) 

^p G L°°(0, T; H\D)), G L°°(0, T; H^D)), (3.24) 

^i/iOp-i g p,i0(f2)), ^Vp-^ G l2((o, T) X D), (3.25) 

u G L2(0,r;L2(D)),piu G l7(0,T) x D), 

VVp G T°°(0, T; L^(D)), VdA^p G L‘^{0, T; l2(Q)). 


(3.26) 
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In particular, we have 


[ p\V^ logpp 

Jn 


dx dt < C, 


which yields 

\Wp\\L‘2{0,T-,H^{n)) + W'^P'^ llL4(0,r;L4(Q)) < C, (3.27) 

where the constant C > 0 is independent of rj, 5. That is, this inequality is still true after 
r/ —)• 0 and <5 —)• 0. 

Thus, we have the same estimates as in Lemma [3. 2 1 at the levels with r] and 5. Thus, we 
deduce the same compactness for (p^, u^) and {ps, U5). Here, we focus on the convergence 
of the terms and 5pVA^p. Here we pass to the limits with respect to p first, and 

then with respect to 5. Here we state the following two lemmas. 

Lemma 3.6. For any p^ defined as in Provosition 17.71 we have 

rT r 


V 


Pjj dx dt 


0 


as r/ —)• 0. 

Proof. By (|3.22l) . we have 

We notice that 


(ln(—))+ dx < C(ro) < 00. 
1 Pv 


y G 


In(-). 

y 


is a convex continuous function. Moreover, Fatou’s Lemma yields 

f (ln(-))_|_ dx < f liminf(ln( —))+dx 
Jn P Jn Pv 


< lim inf 

rj^O 


/(ln(G)) 

Jfi Pv 


dx, 


and hence (ln(i))+ is in L°°(0, T; L^(H)). It allows us to conclude that 


(3.28) 


|{x : \p{t,x) = 0}| = 0 for almost every t, 

where 1^41 denotes the measure of set A. 

By {Pv)t = —— prjdivu^, and thanks to (I3.23p - (|3.27l) . we have 

{pr,)t G L2(0,T;L3(J^)) + l2((0,T;L2(H)). 

This, together with (j3.24p . up to a subsequence and the Aubin-Lions Lemma gives us that 
converges to p in L^(0, T; L^(II)), and hence pr/ ^ p a.e.. 

Thanks to (13.28^ . we deduce 

pp-^° 0 a.e. (3.29) 

By (I3.25P and Poincare’s inequality, we have a uniform bound, with respect to p, of 

pp-io € L^{0,T-,L\n))nL\0,T-L^in)). 

The PP — interpolation inequality gives 


IIW 


-lOi 
r? I 


L3(0,T;Li(O)) “ llL°°(0,r;Ll(O)) 


3 

\\dPv IIl1(0,T;L3(o)) — 
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and hence rjp^ is uniformly bounded in La (0, T; La (n)). This, with (l,'1.29h . yields 

r]p~^^ ^ 0 strongly in L^(0,T; L^(n)). 


□ 


Lemma 3.7. For any ps defined as in Provosition \3.1\. we have, for any test function (p, 

rT 


1 h 


psVps(p dx dt 0 


as (5 —)• 0. 

Proof. By (|3.24l) and (I3.26p . We have uniform bounds with respect to 5 of 
psG L°°{ 0 ,T-L^{n)), V 6 psGL°°{ 0 ,T-H\n)), Vdps G L\ 0 ,T; 

This, with Gagliardo-Nirenberg interpolation inequality, yields 

||V%||ia<C||V>5|lil||p5|li|. 

Thus, we have 

/ ^( [ dt<C sup (||p5||L3(f^))9 [ [ psl"^ dx dt, 

Jo \Jn ) t6(o,T) Jo Jn 

which implies 


(5^|V%| G Lf{0,T-,L^{n)). 


(3.30) 


For the term 


-fl 


in 

we focus on the most difficulty term 

sfi' 


'0 Jn 


psVpsip dx dt =—6 / / A'^div{pi(p)A^PS dx dt, 

Jo Jn 

Cip>) [ [ V 6 \V^^ps\d^\V^p 5 \S^ dxdt 

Jo Jn 


< 


A^iy ps)A^ psp> dx dt 

< C{ip)5^ II v^V^°P5||2,2(0,t;L2(o)) Il<5i9 V®P5|| ^ w 

^ 0 


L‘5'(0,T;L3(r2)) 


as 5 —)• 0, where we used (j3.30p . 

We can apply the same arguments to handle the other terms from 

rT 


6 f [ A‘^d\\{pspi)A^psdxdt. 

Jo Jn 


Thus we have 


as d —^ 0. 


/o Jn 


psVA^psif dxdt ^ 0 


□ 


Here we have to remark that p.27p is still true even after vanishing p and 5. Thus, 
letting r/ —)• 0 and 5 —>■ 0, we have shown that {p, u) solves dni). 

Meanwhile, due to weak lower semicontinuity of convex functions, we have do]) by 
vanishing p and 6 in energy inequality ()3.2ip . Similarly, we can obtain BD-entropy (II.4h 
by passing into the limits in (I3.22p as ?? —)• 0 and <5 —)• 0. 
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3.4. Other Properties. The time evolution of the integral averages 


te{o,T) 


In 


{pu)(t, x) ■ 'ip{x) dx 


is defined by 
d 
dt 

+ 


/ {pu.){t,x) ■'ip{x) dx = / /9U (g) u : dx + / p'^dw'ipdx+ / pDuS/'ipdx 

Jo. Jo. jQ Jo. 

ro / u'ljjdx + n / p\u\‘^u'ijj dx + 2 k / ^/p'ljj dx + k / Ay/p^/pdivijj dx. 

J ^ J ^ J ^ J ^ 


(3.31) 


All estimates from (m and (11.41) imply (|3.31l) is continuous function with respect to 
t G [0,r]. On the other hand, we have 


and hence 

We notice 

thus 


pu G L°°(0, T; L2 (11) n L^(0, T; L'^in)), 


pneCi[0,T]-,Ll^^,m. 


{Vp)t = “2 • u, 


ll(^/p)^llL2((o,^)xn) < C'||-v/pdivu||2,2((o,T)xO) + llL4((o,r)xQ)IIP'‘u||i4((o^'r)xO)- 

This, with Vy/p^ G L°°(0, T; L^(ll)), we have 

y/p strongly in L^(0,T; L^(n)). (3.32) 

Because 

{pu.)t = —div(/9U (g) u) — Vp'*' + div(/?Du) — xqu — npluPu + KpV (—^ ) , 

\ Vp J 

thus we have {pu)t is bounded in L^(0, T; 1T“^’^(0)). Meanwhile, we have 

V(pu) = (pLu) • Vy/pp^ + y/p^Vu, 

which yields V(pu) G L^(0, T; Ls (H)) + L^(0, T; L2 (H)). The Aubin-Lions Lemma gives 
us 

PkUk ^ pu strongly in L^((0,T) X n). (3.33) 

Applying Lemma 13.31 with (I3.32|) , (I3.33|) , and 

rT 


[ [ P.K 

Jo Jn 


' dxdt < C < 00, 


we have 


y/^u^ y/pu strongly in L^(0, T; L^(0)). 
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